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ABSTRACT 

We consider the prospects for measuring the pairwise kinetic Sunyaev-Zel'dovich (kSZ) signal from 
galaxy clusters discovered in large photometric surveys such as the Dark Energy Survey (DES). We 
project that the DES cluster sample will, in conjunction with existing mm-wave data from the South 
Pole Telescope (SPT), yield a detection of the pairwise kSZ signal at the 8-13tr level, with sensitivity 
peaking for clusters separated by ^100 Mpc distances. A next-generation version of SPT would allow 
for a 18-30cr detection and would be limited by variance from the kSZ signal itself and residual thermal 
Sunyaev-Zel'dovich (tSZ) signal. Throughout our analysis we assume photometric redshift errors, 
which wash out the signal for clusters separated by <50 Mpc; a spectroscopic survey of the DES sample 
would recover this signal and allow for a 26-43er detection, and would again be limited by kSZ/tSZ 
variance. Assuming a standard model of structure formation, these high-precision measurements 
of the pairwise kSZ signal will yield detailed information on the gas content of the galaxy clusters. 
Alternatively, if the gas can be sufficiently characterized by other means {e.g. using tSZ, X-ray, or weak 
lensing) , then the relative velocities of the galaxy clusters can be isolated, thereby providing a precision 
measurement of gravity on 100 Mpc scales. We briefly consider the utility of these measurements for 
constraining theories of modified gravity. 



1. INTRODUCTION 

Measurements of the peculiar velocities of massive 
halos, in conjunction with expansion history measure- 
ments, can be used to constrain cos mological models and 
test gravity on cosmological scal es ( |Kosowsky fe Bhat- 
|tacharya||2009| Reyes et al.||2010 ). The kinetic Sunya ev- 
Zerdovich (kSZ) effect ( |Sunyaev fc Zerdovich||1972j ) ef- 
fect provides a means to measure peculiar velocities. 
CMB photons gain or lose energy after Thomson scatter- 
ing on free electrons moving with respect to the CMB, 
and the resulting shift in the CMB temperature in the 
direction of the electrons is known as the kSZ effect. 
The kSZ effect w as recently measured at the 3.8cr level 



2012) using mm-wave data from the Ata- 
ogy Telescope (A CT) in conjunction with 



( Hand et al. 
cama Cosmo' 

the BOSS spectroscopic catalog (SDSS-III Collaboration 
et al.|2 012). This was a measurement of the pairwise kSZ 
signal; two massive halos will tend to fall towards each 
other, and the resulting CMB temperature difference is 
the pairwise kSZ signal. 
The cosmological utility of a catal og of kSZ-derived 



velocities was s tudi ed in detail in Bhattacharya fc 



Kosowskyj ( |2007[ ) and |Kosowsky fc Bhattacharya[ (2009 ) 



and the prospects for measuring growth using the angular 
correlation of galaxies found in the Dark E nergy Survejj^] 
(DES) was studied by |Ross et~aL] ( |2011| . In this letter 
we consider the prospects for detecting the pairwise kSZ 
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signal using existing or fut ure mm-wave da t a from the 
South Pole Telescope (SPT, |CarIstrom et al.| ( [2011] )) and 
galaxy clusters discovered by the DES, which is coming 
online at the time of writing. We simulate photometric 
redshift errors and realistic mm-wave data, including the 
noise introduced by the tSZ and kSZ effects, and find that 
these datasets will yield strong detections of the pairwise 
kSZ signal. 

2. ANALYSIS 



Our analysis re lies on the simulations presented in Se- 
hgal et al. (2010), hereafter S10. This work constructed 
simulated maps of the mm-wave sky with contributions 
from the CMB, the thermal and kinetic SZ effects, and 
emissive sources. The pieces of S10 that are relevant for 
our analysis are the halo catalog and the simulated maps 
of the kSZ and tSZ signals. The SZ signals were cal- 
culated per cluster in post-processing; the gravitational 
potential was calculated according to the N-body parti- 
cle density, and the gas density and pressure were calcu- 
lated assuming hydrostatic equilibrium and a polytropic 
equation of state. S10 calculated the kSZ signal after 
assigning a single, mass-averaged velocity to all gas as- 
sociated with each cluster. We scale the S10 tSZ map 
by a factor of 0.7, which is in broad agreement with the 
signals me asured in recent tSZ cluster surveys (Benson 
et al. 201 ip and m easurements of the tSZ power spectrum 
i Reichardt et al.||2012 |. 



2.1. Cluster Sample 

We construct two DES-like cluster samples from the 
S10 halo catalog. The first, dubbed DES-frase, includes 
clusters with 1 x 10 14 M m < M 5nnr < 3 x 10 14 M n and 



0.2 < z < 1.0. The lower mass threshold, M 50 o c 
1 x 10 14 Mo, correspo nds to approximately 20 galax- 



ies (Rykoff et al. 2012), which is the expected selection 
threshold tor DES clusters. We ignore scatter in mass (at 



2 



fixed optical richness) and offsets between the optically- 
derived cluster center and the center of gas mass, both 
of which would need to be modeled in a more realistic 
analysis. The second, more aggressive sample, dubbed 
DES-agg, has the same redshift range but a lower mass 
range: 5 x fO 13 M Q < M 500c < 1.5 x 10 14 M . Such low 
mass systems might be obtained in exchange for lower 
purity and higher mass scatter. 

We have limited the upper mass range of these samples 
to create more uniform samples and to mitigate strong 
tSZ signals. This requirement reduces the number den- 
sity of clusters by ~10% per sample, leading to ~4 (16) 
clusters per deg 2 in the DES- base (DES-agg) sample. 

We assign photometric redshift errors to the cluster 
redshifts according to a z = 0.01 x exp(z/2.5), approx- 
imating the DES photometric errors presented in The 
Supplements for the Dark Energy Survey White Papei^ 
The mean redshifts are the redshifts provided in the S10 
catalog, which do not include distortions due to velocities 

2.2. mm-wave Data 

We simulate observations of the mm-wave sky by the 
SPT. We consider two surveys: 

• SPT-SZ - the initial, 2500 deg 2 survey, for which 
we assume white noise levels of [40, 18, 80] /xK- 
arcminute and beam widths of [1.6, 1.05, 0.85] ar- 
cminute FWHM at [95, 150, 220] GHz, and 

• SPT-3G - a future, third-generation survey, which 
covers the same 2500 deg 2 area, and for which 
we assume white noise levels of [4.3, 2.5, 4.3] /iK- 
arcminute and beam widths of [1.6, 1.05, 0.85] ar- 
cminute FWHM at [95, 150, 220] GHz. 

While each of these surveys has or would have cov- 
erage at 95, 150, and 220 GHz, we use only the 150 
and 220 GHz data in the simplified analysis pres ented 
here (with one exception, discussed in Section 2.5). We 
assume a background of randomly-distributed, emissive 
sources with y^C^ = 7 (21.5) /iK-arcminute at 150 (220) 
GHz that is perfectly spatially correlated be tween 150 
and 2 2 GH z, consistent with the results of Rcichardt] 
et al.| ( |2012[ ). We optimally combine the 150 and 220 
GHz data as a function of multipole to minimize the to- 
tal power due to instrumental noise and emissive sources, 
and the resulting noise power is defined as Ng. The 
beam-deconvolved noise p ower a t £ — 4000, where the 
kSZ sensitivity peaks, is \J -/V4000 ^22 /iK-arcminute (6 
/iK-arcminute) for SPT-SZ (SPT-3G). 

We simulate mm-wave data by combining the kSZ and 
tSZ skies from the S10 simulations with Gaussian real- 
izations of Ni, and Gaussian realizations of Cf MB , the 
best-fit ACDM CMB power sp ectrum as constrained by 
SPT+WMAP in |Keisler et aT[ ( |201l| ). We use n sido = 2 13 
Healpis]^] maps throughout our analysis. 



2.3. Signal Estimation 

https:/ /des. fnal.gov/survey_documents/DES- 
DETF/Supplements_DES-DETF_vl . 6.pdf 
8 http://healpix.jpl.nasa.gov 



Next we estimate the pairwise kSZ signal of the DES- 
like cluster samples in the simul ated mm-wave data . One 
possible technique, adopted by Hand et al. (2012), is to 



convolve the mm-wave map with a roughly optimal spa- 
tial filter, and, for all cluster pairs {ij} within some range 
of pair separation, measure 



Pksz(r) 



^i<j(Ti — Tj)ci 



(1) 



where Ti is the filtered temperature at the location of 
cluster 2, Yi is the comoving distance to cluster i, 

r — Vj , and Cij = f ,y • r '"^ r3 is a factor to account for the 
on-sky projection of the clusters' separation. 

We instead adopt a harmonic-space technique, which 
we suspect may be equivalent to the real-space technique 
described above. We assume that, for all pairs of clusters 
{ij} within a bin of cluster separation, the gravitational 
attraction of cluster j induces, on average, a kSZ signal 
at the location of cluster i given by 



(2) 



T l (6)=T c ij 1 



where 9 is the angle from cluster i and 9 C is the an- 
gle subtended by r c =0.36 Mpc (physical distance). In 
other words, we assume a /3-modcl profile for the optical 
depth to Thomson scattering with (3 = 1 and r c =0.36 
Mpc, which approximates the optical depth profiles of 
the DES-like clusters in the S10 simulations. We account 
for projection effects using the c^- factors, which we cal- 
culate using the clusters' angular positions and simulated 
photometric redshifts. 

Our technique boils down to fitting for T$, the over- 
all normalization, as a function of cluster separation and 
redshift. In practice this amounts to populating a blank 
Healpix map with the projection factors at the loca- 
tions of all clusters in a bin of cluster separation. Specif- 
ically, we populate the position of cluster i (j) with 
(-Cij), and count each pair exactly once. We convolve 
this map by the /3-model profile to yield the kSZ tem- 
plate map, )C, and its spherical harmonic components 
ICi m . Finally, we estimate Tq as 



Re[d em K.* em ]M; 



(3) 



where d^ m is the spherical harmonic component of a 
simulated mm-wave map and M e = N e + C^ MB accounts 
for all sources of noise power. We repeat this procedure 
for several bins in pair-averaged redshift, z avg , and clus- 
ter pair separation, r com . There are three redshift bins, 
[{0.2, 0.4}, {0.4, 0.7}, {0.7, 1.0}], and five pair separa- 
tion bins, [{10, 30}, {30, 70}, {70, 130}, {130, 200}, {200, 
300}] Mpc. The kSZ signal in the bin of largest separa- 
tion may be affected by the 1000 comoving Mpc/h size 
of the simulation box used in S10. 

Many of our assumptions (e.g. that all clusters within 
a given range of redshift and cluster separation have an 
identical gas profile) are surely inaccurate at some level, 
but these assumptions can be regarded as definitions of 
the analysis technique. The technique is then calibrated 
by running the same analysis on simulated skies, as was 
done here. The challenge, of course, will be modeling 
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and marginalizing over the inadequacies of the simulated 
skies. 

2.4. Noise Estimation 

The S10 SZ maps are duplicated per octant of sky, 
and we assign one 2500 deg 2 SPT survey area, includ- 
ing unique photometric redshift errors, per octant. The 
CMB, instrumental noise, and point source background 
are generated as Gaussian realizations of Me. The kSZ 
amplitude To is estimated once per sphere, and the scat- 
ter per survey is v8 larger than the scatter per sphere. 

We simulate >120 spheres per data scenario. This en- 
semble of T estimates is used to construct a covariance 
matrix C in bins of cluster separation. C is calculated 
per redshift bin, per data scenario. 

The noise introduced by photometric redshift errors, 
instrumental noise, the background of point sources, and 
the CMB f\ are accounted for naturally in this scheme, 
but there is additional noise introduced by the kSZ and 
tSZ signals. Both of these arise from the finite number 
of clusters in the sample. We emphasize that the tSZ 
signal is not a source of bias; the pairwise kSZ signal is 
differential, and the tSZ will, on average, difference out. 

Ideally the kSZ and tSZ noise would be estimated in 
the same manner as the other sources of noise — per 
octant, per Healpix realization — but we have only one 
independent octant of S10 SZ maps. Instead, we esti- 
mate the kSZ and tSZ noise using a bootstrap technique. 
We divide the cluster sample into eight redshift bins and 
measure the scatter in Tq (at fixed cluster separation) 
across these redshift bins. We translate this result to 
our main analysis by accounting for the different redshift 
bin widths and by assuming that the kSZ and tSZ co- 
variance matrices have the same dimensionless bin-bin 
correlation structure as the main covariance matrix C 
described above. 

One source of noise that is absent in our analysis is 
velocity substructure within a cluster. S10 assigned a 
single peculiar velocity to all gas associated with a par- 
ticular cluster whe n calculating the kS Z map. Howeve r, 



( [20031 ), |Holder| ( [2004] )) 



several works {e.g. Nagai et al 

have shown that velocity substructure limits tne accu- 
racy with which mm-wave data can estimate a cluster's 
velocity. The systematic floor is ~50-100 km/s per clus- 
ter, or 8T^2 fiK for the clusters considered here. This 
is completely subdominant for SPT-SZ noise levels, and 
unlikely to change our results for SPT-3G by more than 
10%. 

We have also neglected uncertainty in a z , the width 
of the photometric error distribution. The effect of pho- 
tometric errors on the kSZ signal is a strong function 
of cluster pair separation and most strongly affects <50 
Mpc separations, as shown in the bottom panels of Fig- 
ure [l] For this reason we expect the marginalization over 
o~ z to only mildly weaken the utility of the pairwise kSZ 
measurement, especially when a prior on o z is provided 
from the broader DES analysis program. 

2.5. Results 

9 We have neglect ed lensing of the CMB by the clusters, but 
|Hezaveh et al.| j2012fr show that this effect only decreases the CMB 
variance at the cluster locations. 



We use the simulations described in the previous sec- 
tions to determine the detection significance of the pair- 
wise kSZ as a function of pair-averaged redshift. Specifi- 
cally, for each redshift bin we compute \ 2 = d T Cd, where 
d is the mean kSZ signal as a function of pair separation 
in this redshift bin and C is the covar iance matrix de- 
scribed above, and calculate SNR = ^/EiX 2 , where the 
sum is over the three redshift bins. We summarize the 
results in Table [T] and Figure [ll 

We find that existing data from the SPT-SZ survey 
will, when combined with either of the DES samples con- 
sidered here, result in a strong detection of the pairwise 
kSZ signal. The DES- 6ase (DES-agg) sample yields a 
8.2it (12.9<r) detection. The smaller optical depths of the 
clusters in the DES- agg sample are made up for by the 
increased number density. For both samples, the mea- 
surement is largely limited by instrumental noise in the 
SPT-SZ data. 

As expected, we find that data from a future, third- 
generation survey, SPT-3G, would result in even more 
precise measurements of the pairwise kSZ signal. The 
DES-&ase (DES-agg) sample yields a 18<r (30cr) detec- 
tion. In contrast to SPT-SZ, the SPT-3G constraints are 
limited by kSZ and tSZ noise, with roughly equal contri- 
butions from each. 

We considered using a linear combination of 95 and 
150 GHz SPT-3G data to eliminate the tSZ signal and 
its associated variance. We found that the elimination of 
tSZ noise was not worth the effective increase in instru- 
mental noise; the resulting detection significance on the 
DES-agg sample was 18er, considerably weaker than the 
30c obtained using the original analysis. 

Finally, we explored how these constraints would 
change in the event of spectroscopic followup of the DES 
cluster sample, such as would be possible with DESpec 
(Abdalla et al. 2012). We assume a spectroscopic red- 
shift error of a z = 0.002 per cluster and find that SPT- 
SZ (SPT-3G) would yield a 12.5cr (26cr) detection on the 
DES- base sample and a 26c (43c) detection on the DES- 
agg sample. 

2.6. Implications for Testing Gravity 

The pairwise kSZ signal directly measures the average 
pairwise ionized momentum of a given cluster sample. 
The results of the previous section demonstrate that pre- 
cision measurements of this signal should be available in 
the coming years, and these have the power to yield de- 
tailed information on the gas content of the DES clusters. 

To go beyond gas physics — to constrain cosmology 
or gravity — we must convert to average pairwise veloc- 
ity. This requires an independent measurement of the 
average ionized gas profile, say from tSZ or X-ray obser- 
vations. Additionally, if there is a well motivated con- 
nection between the total cluster mass and the mass in 
ionized gas, then measurements of the mean mass of the 
sample through weak lensing or angular clustering can 
also constrain the conversion from kSZ signal to velocity. 
For the remainder of this section we assume that the un- 
certainty in the gas as constrained by independent data 
is subdominant to the uncertainty on the pairwise kSZ 
measurement. 

The mean pairwise velocity of halos (vh)(r) as mea- 
sured through the kSZ effect as function of separation r 
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TABLE 1 
Sensitivity to pairwise kSZ signal 





DES-6ase (photo-z) 


DES-agg (photo-z) 


DES-6a.se (spec-z) 


DES- agg (spec-z) 


SPT-SZ 


8.2cr (12%) 


12.9cr (7.8%) 


12.5(7 (8.0%) 


21(7 (4.8%) 


SPT-3G 


18cr (5.5%) 


30cr (3.3%) 


26cr (3.9%) 


43(7 (2.3%) 



The projected detection significance of the pairwise kSZ measurement for different combina- 
tions of CMB data and cluster samples. 



DES-BASE, photo-z DES-AGG, photo-z 




Comoving Separation (Mpc) Comoving Separation (Mpc) 



DES-BASE, spec-z DES-AGG, spec-z 




Comoving Separation (Mpc) Comoving Separation (Mpc) 



Fig. 1. — The simulated pairwise kSZ signal for the DES-6a.se (left panels) and DES-agg (right panels) samples, and for photometric 
(top panels) and spectroscopic (bottom panels) redshifts. The black curve shows the mean signal, while the error bars show the projected 
uncertainties using data from SPT-SZ (blue, left offset) and SPT-3G (red, right offset). The dashed, gray curve in the bottom panels 
shows the mean signal using photometric redhifts, to highlight the loss of signal at r com < 50 Mpc. The errors include contributions from 
instrumental noise, a background of point sources, the CMB, the kSZ signal itself, and residual tSZ signal, and are mildly correlated between 
bins of cluster separation. The second vertical axis shows the corresponding velocity per cluster, assuming optical depths of r = 0.0025 for 
DES-6a.se and r = 0.0017 for DES-agg, which are typical of the clusters used in each sample. 
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is in linear perturbation theory given by ( [SchmidtpOlot 

(v h )(r)=2b£ Sv (r). (4) 

Here, b is the kSZ signal- weighted mean of (b— 1) of halos 
in the given mass bin, b is the linear Eulerian bias, and 
£sv (r) is the linear cross-correlation between density and 
velocity. If the kSZ signal is simply proportional to gas 
mass which is proportional to halo mass, b reduces to the 
mass- weighted average of b — 1. We assume that there is 
no velocity bias of halos, which holds on the large scales 
of interest. 

The velocity-matter cross correlation is given by 



2ir z 



dk kP(k, z)jx(kr), 



(5) 



where / = d\n D /d\na is the growth rate, P(k,z) is 
the matter power spectrum at redshift z, and ji denote 
spherical Bessel functions. 

Consider a modified gravity (MG) model that leads 
to a scale-independent modification of growth on lin- 
ear scales (one example is the Dva li-Gabadadze-Porrati 
(DGP) model ( |Dvali et al. |2000| and its generaliza- 
tions) . The effect ot such a modification of gravity on 
the velocity-matter cross-correlation at fixed redshift is 
given by 



£<5u,ACDm('", z) 



[fcrl(z)] M G 
L/^IMIacdm' 



(6) 



where cr|(z) is the power spectrum normalization at red- 
shift z. On linear scales, the mean pairwise velocity as 
measured through kSZ thus constrains &/<r|. Indepen- 
dent knowledge of b can then be used to turn the kSZ 
measurement into a constraint on /<t|. When includ- 
ing constraints on the expansion history and primordial 
amplitude of fluctuations, this can then be turned into 
a constraint on modified gravity. Specifically, for the 
normal-branch D GP model with ACDM expansion his- 
tory considered in Schmidt | ( |2009 ) which is parametrized 
by the cross-over scale r c , a measurement of (r) with 
an overall SNR of 10 (20) — as should be possible using 
the DES cluster sample in conjunction with existing data 
from SPT-SZ (future data from SPT-3G) — would yield 
a 95% confidence level lower limit of r c > 3000 (7000) 
Mpc, assuming that the expansion history, halo bias, and 
kSZ- velocity relation are known perfectly (Fig. [2]). This 
is an interesting constraint, since a modified gravity ex- 
planation for the accelerated expansion of the Universe 
suggests r c ~ 3000ft. -1 Mpc as a natural value in these 
models. 

We note that uncertainty in the cluster gas affects only 
the global normalization of the kSZ signal considered 
here, not the shape of the kSZ signal as a function of sep- 



aration. This can be used to test the scale-dependence 
of gravity on sufficiently large scales, independently of 
gas physics. A well-studied class of m odified gravity 
mode ls which includes the f(R) model (Carroll et al 



20041 invokes a massive additional degree ot freedom. 
The modifications to the gravitational force and veloc- 
ities are generically scale-dependent and suppressed on 
scales larger than the Compton length of the field, typ- 
ically 30 Mpc or less (see Fig. [2]). For this reason, only 
cluster samples with spectroscopic redshifts are useful 
for constraining f(R). Our most optimistic scenario, in 
which SPT-3G observes the spectroscopically followed- 
up DES-agg cluster sample, results in a 15a (6.7%) de- 
tection of pairwise kSZ signal in the the 10-30 Mpc sep- 
aration bin. For comparison, the maxim um currently- 
allowed value of the field amplitude in the Hu & S awickil 
(|2007|) f(R) model is \f R0 \ = 10" 4 ( [Schmidt~etlil.||2009| 
|Lombriser et al.|2012 ) , and results in a ~257o increase m 
£g v over this range. The SPT-3G+DES-ag(?(spec-z) sce- 
nario would therefore allow a 4cr test of this hypothesis. 

3. DISCUSSION 

We have estimated the sensitivity of mm-wave data 
from SPT to the pairwise kSZ signal of clusters from the 
DES. We find that existing SPT data will allow for an 
8-13cr detection of the kSZ signal. If the kSZ-velocity 
relation is known perfectly, this corresponds to an 8-12% 
constraint on /erf, which is comparable to the current 
best constraint on growth from redshift space distortions 
(RSD), an 8.2% constraint on fa s from BOSS < |Reid et al. 
2(512} . 

The difference is that the bulk of the RSD sensitiv- 



ity arises from quasi-linear scales ( Reid fc WhitepOll 1, 
whereas the SPT+DES kSZ signal is most sensitive to 
significantly larger, more easily modeled, linear scales. 
This is due to the DES photometric redshifts, and spec- 
troscopic followup would boost the kSZ SNR by a factor 
of 1.5. A third-generation SPT survey would improve 
any of these constraints by a factor of ~2.1, and is lim- 
ited largely by noise from the kSZ and tSZ signals. These 
measurements have the potential to provide interesting 
constraints on theories of modified gravity. 
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